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Abstract 

We demonstrate that to calculate the self-energy of a heavy quark in the heavy 
quark limit (or the heavy fermion limit in what is called the Baryon Chiral Pertur- 
bation Theory), the use of standard dimensional regularization provides only the 
quantum limit: opposite to the heavy quark (or classical) limit that one wishes to 
obtain. We thus devised a standard ultraviolet /infrared regularization procedure 
in calculating the one- and two-loop contributions to the heavy quark self-energy 
in this limit. Then the one-loop result is shown to provide the standard classical 
Coulomb self-energy of a static colour source that is linearly proportional to the 
untraviolet cutoff A. All the two-loop contributions are found to be proportional 
to Aln(A/A) where A is the infrared cutoff. Often only the contribution from the 
bubble (light quarks, gluon and ghost) insertion to the gluon propagator has been 
considered as the 0(qs) correction to the Coulomb energy to this order. Our result 
shows that other contributions are of the same magnitude, thus have to be taken 
into account. 
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1 Introduction 



Recently, many discussions have been devoted to the definition of the of heavy quark mass. 
This quantity is the basic parameter in the heavy quark physics. Generally (not only for 
the heavy quarks in the heavy quark limit in which we are interested, in the following) one 
of the most popular choices within the perturbative scheme is the pole mass. In reference 
m it has been shown that the pole mass is gauge independent and infrared finite at the 
two-loop level. It was checked at the three-loop level in where the relation between 
the mass of the modified minimal subtraction scheme and the scheme-independent pole 
mass was also obtained. Later authors of work IQ argue that precise definition of the 
pole mass in the heavy quark limit may not be given once nonperturbative effects are 
included (a more analytical approach leading to the identical conclusion has been given 
independently in [^]). To demonstrate that, they adopted the standard effective running 
constant (as a typical non-perturbative quantity) obtained from the summation of leading 
logarithms to calculate the Coulomb energy of a static (or infinitely heavy) colour source. 
This has produced an auxiliary singularity for small momentum (the so-called, infrared 
renormalon) which, they claimed, leads to the uncertainty of order Kqcd/^q in the 
definition of the pole mass. An important observation which served as one of the bases 
to this demonstrative calculation is that the widely used dimensional regularization may 
not be used to analize self-energy diagrams in the heavy quark limit, in particular, for 
separating the infrared and ultraviolet domains. The best example, as we shall also 
present in the following section, is that this regularization method leaves out the linear 
divergence that defines the self-energy of a static quark in the lowest order: the Coulomb 
energy. Also it should be emphasized here that the Heavy Quark Limit does make sense 
only after introducing a scale parameter to the theory with which the quark mass may 
be compared, and there is something awkward about it in the dimensional regularization 
method. 

It may be useful to point out here that a quite analogous situation to the heavy quark 
limit can be found in the problem of the classical limit % ^ {] for the self-energy of an 
electron in QED It was found that in the limit ^ 0, or more precisely for the electron 
Compton wave length h/mc <^ tq, where parameter tq regularizes small distances, the 
energy of the Coulomb field by a change of radius tq is reproduced in the second order of 
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perturabative theory whereas the higher order contributions were shown to vanish in this 
hmit. 

Here we first study the quantum and classical limits for the lowest order contributions 
to the quark self-energy and quark-gluon vertex, and show that the standard dimen- 
sional regularization method only reproduces the quantum limit. Since one may clearly 
identify the classical limit to be equal to the currently fashonable heavy quark limit, we 
could conclude that the dimensional regularization method is not suited for studying this 
case. Then in the following section we calculate the two-loop contributions to the quark 
self-energy in the heavy quark limit, exploiting the conventional ultraviolet and infrared 
regularization method. We have found that each contribution to this order is all propor- 
tional to A In (A/ A) where A and A are the ultraviolet and infrared cutoffs, respectively and 
that, unlike in QED, these contributions do not mutually cancel out. We point out that 
the gluonic bubble diagrams (Fig.4d) (implicitly) taken often as the majour contribution 
in some processes is of the same magnitude as those from vertex and quark self-energy 
diagrams (Figs. 4a-4c) for the self-energy in the heavy quark limit. In the last section 
some comments will be given regarding the alleged infrared renormalon and the Baryon 
(or heavy fermion) Chiral Perturbation Theory. 

2 One-loop diagram 

The self-energy of quark with mass m in the second order of perturbative theory is defined 
by the diagram in Fig.l, 

5m, = S,(m) = ^i^^.V(l^(^^';;^3^7.ly^.^ (1) 

where Cp = {N^ — l)/2Ai'c, for a gauge group SU{Nc). First, we note that this expression 
coincides with the electromagnetic self-energy except for the trivial factor CfqI- While 
the integral has no infrared divergence, it develops the linear ultraviolet divergence. We 
are interested in this quantity in the heavy quark limit: m —>■ oo. As was mentioned in 0, 
the dimensional regularization cannot serve this purpose. To make this point more clear, 
we calculate the quark self-energy by way of introducing an ultraviolet regularization. 
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1 F{-e/K') 

Here, the function fc^/A^) is assumed to decrease rapidly in the Euchdean region 
kl^ = —k'^ +00 providing the uhraviolet convergence of the Feynman integrals. Now, 
the heavy quark limit (or classical limit) means the case where m 3> A in comparison 
with the quantum limit: m <^ A. The both limits may be realized using the Mellin 
representation for the form factor F(— /c^/A^) under calculation of the integral in eq. (1). 
We can write in the Euclidean region 

nki) = f:cn{-nk%r = ^ J ^<oi&, (3) 

where c(n) = c„ for n being positive integers, and < 5 < 0.5. 
Finally, we have 



— <5+ioo 



C(m2 - (k+py) 



2 1 V" dC 1 r(2-C) K ( d^k 4 -2(1 -a) 

= ^''^"^21 J d^^^^^A^rXW) r"^^ - y (2.)^z (aW - F)^-C 

-5+ioo ^ V / V / 

3«3^ 1 "V" dc ,,, r(i + 20 ^m^^C r(l- 0(1 + ... 
= ^^^"^2z y ^^^^^^^UJ r(3 + o • 

—S+ioo 

To get the leading term in m/ A in the quantum limit, we move the integration counter 
to the right and take into account of the first double pole at C = 0. Assuming that 
F(0) = 1, we have 

c qu 3a, „ r, A 
orrio 



C^mjln— + 0(1)1. (5) 



47r 

To get the leading term in A/m in the classical limit, we move the integration counter 
to the left and take into account the first pole at C = —0.5. From the standard relation 
between the Mellin and its inverse transformations 
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we have 



oo 

c(-0.5) =2 JdtF(t^), 





oo 

5m'^ = ^CMc(-0.5) = ^CpA [dtF(t^). (6) 



If we compare the above result with the calculation of the expression in Eq.(l) devising 
the dimensional regularization, 

where A^^^ = 1/e + ln47r — 7, (7: the Euler-Mascheroni constant) and /x is the subtrac- 
tion point, one can see that it only reproduces the quantum limit Eq.(5) for the quark 
self-energy by identifying /i as the ultraviolet cutoff. Thus, the naive dimensional regu- 
larization does not allow to see the heavy quark limit. 

Ordinarily, studies of heavy quark limit may be carried out by using the expansion of 
quark propagator in the inverse quark mass: 



[1 + i/2m + /(A;Vm2)][fc4 + 2rmf{kym'^)] 
iMo _ , 1 



- + 0{-) (8) 



k4 + ie m 



Here, = -iko, Mq = (1/2) (7 + 7°), and f{x) = {l/2){^/T^ - 1). It is readily seen 
that the Dirac spinor u(p) for quark may be chosen as 

' I \ 



u(p} = 







so that u{p)u{p)—u(p)j^u{p)—l, and u(p)^u{p)—0. 

In the second order: Eq.(l), upon introducing the ultraviolet regulator to the otherwise 
linearly divergent expression, the leading contribution reads 

Sm^^ = S2(m) = 27ra,CM /t^S^ (9) 



which obviously coincides with the Eq.(6). One can see that this expression is the classical 
Coulomb energy of a static color source. To give more physical meaning to Eq.(9), one 
may introduce the charge density defined as 



Pin 



(2^ 



F{k'^) 



(10) 



supposing that F(P) > 0, and the electron (quark) "radius" Tq = 1/A. 
Then, we have 
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47r Iri — 



P{ri) 



(11) 



in accordance with the classical result. 

As was mentioned above, there is no infrared divergence in the second order, but it 
plays the crucial role in higher orders in perturbative expansion. 

To understand its role, we compute the quark-gluon vertex induced by both quark- 
gluon and self-gluon interactions: Figs. 2-3. For simplicity we omit all the common coef- 
ficients and color factors. 

(a) Quark-gluon interaction (diagram in Fig. 2): 

By introducing the gluon mass A serving to regularize the infrared singularity we get 



Ar(P,P) 



d^k 1 
{2ti)H A2 - k^ 



— {k + pY 

Using the form factor as in Eq. (2) one obtains 



(12) 



Ai(p,p) = 



21n^-lnS + 0(l)| ifA>m 



1 



7 (4^)2 ^ ^ 



2 In f + 0(11 



if A < m 



One can see that the infrared singularity defines the behavior of the vertex in the heavy 
quark limit. Again result obtained from the dimensional regularization for ultraviolet 
divergence corresponds only to the quantum case: 



m 



AT, + In — + 0(1 
//2 



(13) 
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Note that if we do not introduce the infrared cutoff within the dimensional regularization 
scheme, we find the auxihary 1/e pole in the expression for the vertex in Eq.(12) leading 
to 

^''^M = (4^{-2iV/^ - Nr + 31n ^ + 0(1)}. 

(b) Self-gluon interaction (diagram in Fig. 3): 

This contribution may be written as 



This integral turns out to be free from infrared singularity. Using the ultraviolet form 
factor gives 



' rw{31nS + 0(l)} ifA»m 



A2(P,P) = < 

.0(A) ifA«m 

It is readily seen that this diagram gives no contribution in the heavy quark limit. Again 
result obtained from the dimensional regularization for ultraviolet divergence corresponds 
only to the quantum case: 

A^(P,P) = ^'^(^{SA^^ " ^^"^ + ^^^^}- ^^^^ 

3 Two-loop Diagrams in the Heavy Quark Limit 

The contributions to the quark self-energy from the two-loops (or fourth order in coupling 
constant) are defined by the diagrams in Fig. 4. As a consequence of our discussion in 

the previous section, we shall not adopt the dimensional regularization, but merely use 
the ultraviolet form factor Eq.(2) together with the infrared regularization of the gluon 
propagator in the following. 
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(a) Diagram in Fig. 4a. 

With A and A being the ultraviolet and infrared cutoff masses, as in the previous 
section, one finds 



where = A^A^; C, 
I {k^) is equal to 



(a) 



6m'f^ = Ei'^V) = ^^.'A#(«') (16) 
1 (for QED), and ftHH'' = -Cf/2Nc (for QCD). The integral 



d^h fd^k2 F{kf) F{kl) 



(27r)V (27r)4 kf + k'^ + {ku + ^e) {ku + A;24 + ie) {k-^A + ie) 

{E) {E) 



(17) 



Here, to be definite, the integration is over the dimensionless Euclidean variables denoted 
as {E) under the integration sign. The calculation of this integral is given in Appendix. 
We have 



(27r)3 ^2 



(18) 



(b) Diagram in Fig. 4b. 

Here one obtains 



5m? = EP(m) = Cf^,^A7f («^) 



(19) 



where C, 



if tH = {Nc/2)Cf (of course there is no QED counterpart to this 



diagram). The integral is given as 



(b) 2 ^ fd^ fd^ r d^ks F{kl) F{kl) F{{ki-k2r) 
^ ^'^ ^ J (27r) V (27r) V (27r)^ kj + kl + «2 (A;i - k2f + 



{ku + ie){k24 + i^) 
with the numerator equal to 



(20) 



AT, = (^1 - 2 }i2)MoYMoY + 7'^o(^i+ h)MoY + l^M^^M^ih -2h)- 

7 



Since the integration is over Euclidean variables here, = 'j^ik^ — 'jk. With this it is 
easy to check that u{p)Nhu{p) = 0, which means that the contribution from the diagram 
in the Fig. 4b behaves as A/m for large m, thus vanishes in the heavy quark limit. 

(c) Diagram in Fig. 4c. 

This contribution may be written as 



6m 



(m) 



(c)„4 A Ac)/ 2 



K 



(21) 



where C^'^ = 1 (for QED), and tHHH^ = Cj, (for QCD). The integral ifin,^) is equal to 



d^ki fd^k2 F{kf) F{kl] 



{2nYJ {2nY kf + k^ + { {k^ + ie) {ku + k24 + ie) {ku 



It 



1 



^4 (22) 
[ku + it) {k2A + it) {kii + it) \ ' 

Here, we have taken into account the mass renormalization in the second order (the dark 
disk in Fig. 4c). From the above result, it is easy to see that in QED the contributions 
coming from Fig. 4a and Fig. 4c cancel. This is a crucial point in solving the problem 
of classical limit — > in QED since the contribution to the self-energy of a massive 
fermion from the photon vacuum polarization is proportional to 



W^{k) oc {g^^^e -k^k") In 

and vanishes in the infinite mass limit: m — oo. 
For QCD one finds 



1 + 



5mP = 27ra,C^A J- 



(Pk F{k 



2\ r 



a.. 



-— 2C^ln(l/K2) 
An 



(23) 



(27r)3 p 

(d) Diagram in Fig. 4d. 

The contribution from the diagrams in Fig. 4d, involving the vacuum polarization by 
massless quarks, gluons, and ghost is written as 



<5mf = Sf (m) = 2ira,C,A [^^^U-{-k') 



(24) 
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In the above expression, the renormahzed vacuum polarization reads (see Appendix 
for a somewhat detailed calculation) 



1 



a 



— I da In 

An 





1 + ail - a) — 



2Nr. 



l + 4a(l 



a 



4n/Q;(l 



a 



Air 



0(1) 



(25) 



upon expanding in small parameter k^. Here the appearance of the coefficient 
b = {lONc — 2nf)/3 is due to the use of the Feynman gauge. 

Summing up all the contributions from the one- and two-loops, one finds 



6m 



27ra^CirA 



J 



d^k F{P) 
(27r)3 P 

d^k F{P) 



All 



NM{l/te) + —\Yi{P/K 



(26) 



l-^ln(AVA2)Ar, 



2 h_ 

3 ^ K 



(27r)3 p 

Here, we express the final result in term of the first coefficient in the Gell-Mann-Low 
function h = {llNc — 2nf)/3. By adopting the number (fiavor) of light quarks to be three: 
Hf = 3, we can see for A^^^ = 3 that the 0{as) corrections to the leading static Coulomb 
self-energy coming from the bubble insertions to the gluon propagator is of the same order 
of magnitude as those from other diagrams: Figs.4a-c, so the latter contribution should 
be included. 

In ref. the self-energy of the heavy quark in the heavy quark limit was calculated 
by introducing the running constant as(^^) in the Coulomb energy expression: Eq.(9), 
and by introducing the explicit integration cutoff in place of the ultraviolet form factor 
F{-e/A^) (Fig.5), 



6m 



Stt j- d^k asij?^ 



|fc|<MO 



(27r)3 p 



(27) 



with 



as{k^) 



{aM)h/Alr)\n{^il/k^ 



f28) 
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According to Ref. [3], the zero of the denominator in the above expression defines the 
position of the infrared renormalon, leading to the uncertainty in the pole mass. Here, 
we are interested in the O^aK/iD) contribution from Eq.(27), 



Here one might regard fiQ as corresponding to our A (in this respect one could argue 
that the value of the integration cutoff in Eq.(27) should be different from the scale /iq 
introduced in Eq.(28). This would introduce an additional contribution of a logarithmic 
type to Eq.(29), which would however not prohibit the appearance of the alleged renor- 
malon). This result may be compared with the part of the expression in Eq.(26) which 
is proportional to b/Nc- Since we have not calculated the higher loop contributions, we 
cannot draw any definite conclusion. However, in view of the difference in our result: 
Eq.(26) and that of ref. [3], it might be possible that due to the contributions other than 
the bubble insertions in the gluon propagator the characteristics (position, etc.) of the 
alleged infrared renormalon in the heavy quark self-energy would be different from the 
original one, or even the renormalon might not be there. 

4 Discussion 

To summarize our present note, first we have demonstrated that the popular dimensional 
regularization method yeilds only the quantum limit, thus unsuited for the calculation of 
the self-energy of the heavy quark in the heavy quark (or the classical) limit which is just 
the opposite of the quantum limit. Thus we employed next the conventional method for 
ultraviolet and infrared regularizations to calculate the one- and two-loop contributions 
to the self-energy of the heavy quark in the heavy quark limit. The one-loop result is 
found as the standard Coulomb-like self-energy of a static colour souce whose radius is 
characterized by the ultraviolet regularization mass A of the theory. Then all the two-loop 
contributions have turned out to be proportional to Aln(A/A) with the infrared cutoff 
A, while the corresponding contributions in QED internally cancel out to vanish |^. Of 
those two-loop contributions, only the bubble insertion to the gluon propagator in the 
lowrst order (Coulomb) contribution (Fig.4d) has been discussed often. We have found 
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that other contributions are of the same size and have to be taken into account. 

From our study we have reached a couple of observations (possiblly, speculations). 
The first is that since the contributions depicted in Figs.4a-c are not very different in 
magnitude from the bubble insertions to the gluon propagator at the two-loop level, the 
infrared renormalon in the heavy quark self energy resulting from the simple bubble chain 
summation, due to the adoption of the running coupling constant, might be modified 
regarding its position, etc. or might even be absent once a due consideration of the 
remaining higher-loop effects are incorporated. 

The second observation is concerned with the dimensional regularization. Our asser- 
tion regarding the inappropriateness of in this method in the context of the heavy quark 
limit should be valid also in what is called the heavy baryon chiral perturbation theory 
(HBCHPT), (see for example a review and references therein) which combines the 
ordinary Chiral Perturbation theory (CHPT) for the octet peudoscalar mesons coupled 
with baryonic (nucleons, deltas, etc.) degrees of freedom. Here very often the Lagrangian 
(or the S-matrix) for the system is expanded in the inverse power of the baryon mass(es), 
exactly like in the heavy quark expansion realized in the heavy quark limit. The the me- 
son loop integrals arising from the meson-baryon interactions are conventionally carried 
out using the dimensional regularization. Now according to our study in this note, this 
method is only compatible with the opposite quantum limit ! Thus it is our feeling that 
the existing results in HBCHPT should be examined again carefully from the present 
context. 
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Appendix 

(a) The calculation of the integral I^'^\k^). 
Using the obvious equality in Eq. (17) 



ki4 + hi + i€ hi + 
and the well-known formula 



1 

ki4 



ki + ie ki 
one finds 



1 P X 



(27r)3 
where 



We define the function /(t^) as 



{2nf t^ + k^ + 

Using the Mellin representation: Eq.(3) for the form factor F{k'^) and integrating over k, 
one can get 



-5+ioo ^ ''^ 

1 1 rfC 2^c^(l + C)^(l-C) 1 , 1 



— (5+ioo 

for small Finally, we have 
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(a) The calculation of W^''{q'^). 

As an example of the calculation of diagrams in Fig.4, we here present the details of 
the calculation of the vacuum polarization W^^i^q"^) by massless gluons, ghosts, and quarks 
(see, insertions to the gluon propagator in Fig.4d). We proceed in adopting the Feynman 
gauge. We shall use the dimensional regularization in the intermediate calculations, but 
the final result will not depend on the choice of the regularization. 

We define the S-matrix element describing the vacuum polarization as 

S^IJJ dx,dx2Al\{x,)Ul\'^J^{x^ - X2)All{x2) 
with its Fourier transform defined as 

= J dxe-^^^UZlix). 

For simplicity we do not introduce new notation for this value. 

The contribution coming from the gluons and ghosts may be written in the form 

KllliQ) = glr'^^r'^^ / da /.^^ / ,,Mr J I ' T2]2 

L Jgiuon+ghost ^ J [^-Kj^iy-k^ - a[l- a)q^Y 

where 

A^g?uol = ^--(i?-l){^A;^ + [l-4a(l-«)]g^} 

6-L' + (4L'-6)a(l-a)l, 



ghost 



= g^'^'S^-^k^ + 2a{l - - 2a(l - - q^'q"']. 



First, we show that the apparent non-gauge invariant term proportional to g"^""^ does 
vanish. For this objective we use the identity: 

ae + bq^ 1 f ,^Dapa+{2-D)h 



J l-k"^ - aBa^? 2a8 J 



d^'k- 
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which may be proved by using the following expression, 



d 



1 



1 



and integrating it by parts. Then the coefficient of the apparently non-gauge invariant 
contribution becomes 



da 



d^k 



a{l — a) J [— — a{l — a)q^] 



d^k- 



1 



-l + 4L_|a(l-a) 



r2(r>/2-i) {D-\)v\Di2)\ 



Finally one finds the gauge invariant form for the vacuum polarization including both 
gluons and ghosts 



J gluon+ghost 

with Ilgiuon+ghost being equal to 



ngiuon+ghost(?^) = ^ da/i^" J 



d^k [D - 6 - A{D - 2)a{l - a) 



Note that this expression has both infrared and ultraviolet logarithmic divergences 
in four dimension. To regularize the infrared divergence we introduce the gluon (ghost) 
mass A to the denominator. Then the renormalized vacuum polarization is written as 



ngluon+ghost(?^) ~ ngiuon+ghost(?^) ngiuon+ghost (0) 



i 2 

/ (ialn[l-a(l-a)^]|27Vc[l + 4a(l-a)] 
At: J I 



2 

47r J — q 

4A2 



2 Tn^ — 2 



10 + 8- 
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The calculation of the vacuum polarization by quarks turns out to go along the similar 
line with the result identical to the above one up to a constant factor. The combined result 
is the one in Eq.(25). Note that in this equation the momentum is measured in units of 
A. 
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Fig. 1. The lowest order (one loop) quark self-energy diagram. 

Fig. 2. An 0{as) correction to the quark-gluon vertex by quark-gluon interaction. 

Fig. 3. an 0{as) correction to the quark-gluon vertex by gluon-gluon interaction. 

Fig. 4. Quark self-energy diagrams, to two loops. The mass counter term from the lowest 

order contribution is indicated as a blob in Diagram c. In Diagram d inserted in the gluon 

propagator are the gluon, ghost and light quark loops, respectively. 

Fig. 5. Bubble chain in the gluon propagator causing the infrared renormalon, according 

to Ref. 3. 
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